Circuit model for spin-bottleneck resistance in magnetic-tunnel-junction devices 
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Spin-bottlenecks are created in magnetic-tunnel-junction 
devices by spatial inhomogeneity in the relative resistances 
for up and down spins. We propose a simple electrical circuit 
model for these devices which incorporates spin-bottleneck ef- 
fects and can be used to calculate their overall resistance and 
magnetoresistance. The model permits a simple understand- 
ing of the dependence of device magnetoresistance on spin 
diffusion lengths, tunneling magnetoresistance, and majority 
and minority spin resistivities in the ferromagnetic electrodes. 
The circuit model is in a good quantitative agreement with 
detailed transport calculations. 
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Magnetic tunnel junction (MTJ) devices consist of two 
ferromagnetic electrodes separated by an insulating tun- 
nel junction. The resistance of a device increases when 
the ordered moments of the ferromagnets have antiparal- 
lel orientations on opposite sides of the insulating layer. 
Since the relative moment orientations can be altered 
by a small external magnetic field, MTJ devices can ex- 
hibit strong magnetoresistance. Recent progress in 
the fabrication of MTJ devices with large reproducible 
magnetoresistances, has led to interest in their possi- 
ble utility as digital information storage devices for non- 
volatile random access memories or as read heads. Exper- 
imenters have succeessfully fabricated samples with 
magnetoresistances close to theoretical limits proposed 
by Julliere (JJ^]. The Julliere model ignores the spin- 
bottleneck effects which play an important role in studies 
of giant magnetoresistance in magnetic-metal multilayers 
H and in closely related spin-injection experiments [j^J. 
In MTJ devices for which the resistance of the ferromag- 
netic electrodes is much smaller than the resistance of the 
tunneling barrier, the spin-bottleneck effect indeed plays 
a minor role. For modern random access memory appli- 
cations, however, it is important to make MTJ devices 
with much smaller overall resistance to enable operation 
at fast read access times M. The desirability of low- 
resistance microstructured tunnel junctions ||] raises the 
issue which we address here. 

In this paper we propose a simple electrical-circuit 
model for spin-bottleneck contributions to the resistance 
and magnetoresistance of MTJ devices. We use the 
model to evaluate the resistance of the device, for both 
parallel (R p ) and antiparallel (R A ) ordered moment ori- 
entations. An important figure of merit for MTJ devices 
is the magnetoresistance defined by 
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We show that, while both R p and R increase due to 
the spin-bottleneck effect, the magnetoresistance may in- 
crease or decrease. We have established the accuracy of 
the circuit model, by comparing its predictions with those 
resulting from solutions of the coupled spin-up and spin- 
down transport equations for the inhomogeneous device. 
We find that the circuit model gives remarkably good es- 
timates for the magnetoresistance of the device, even at 
a quantitative level. 

Our circuit model starts from the recognition Jh} that 
the two spin-orientations in the ferromagnetic electrodes 
provide two channels which carry current through the 
system in parallel. If the number of up and down spin 
electrons were separately conserved throughout the de- 
vice, it would be described by a circuit model with up 
and down spin channels in parallel and three resistors in 
series for each channel as illustrated in the top panels 
of Fig. |l| and ^. To simplify notation we assume here 
that the ferromagnetic electrodes are identical in length 
(If) and crossection (A) and made from the same ma- 
terial. When ordered moments are parallel, the up-spin 
channel will have the majority spin electrode resistance, 
and the down-spin channel will have the minority spin 
electrode resistance, on both sides of the tunnel junc- 
tion. For antiparallel ordered moment orientations, we 
adopt a convention where the up-spins are in the ma- 
jority on the left side of the junction and in the mi- 
nority on right side. In the Julliere model |^,^], which 
appears to be accurate for most materials of interest, 
the tunnel-junction conductance for each spin channel is 
proportional to the product of factors, one for each elec- 
trode, which are dependent on both the density-of-states 
and the character of the wavefunctions of that spin at 
the Fermi energy. Writing these factors symbolically as 
t a N a and UNi for majority and minority spins, respec- 
tively, we see that when the ordered moment orientations 
are antiparallel the tunneling resistances in the two-spin 



channels are equal: R^ = R 
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For parallel magnetizations Rj,' = R^ a — (t a N a )~ 2 and 
R p '^ = R% = (tiNi)~ 2 . (In these equations N a and N t 
are the majority and minority spin densities of states and 
the factors t a and ij account for the dependence of tun- 
neling amplitudes on wavefunction character.) Julliere's 
formula for the magnetoresistance, which in our notation 
takes the form 
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results from the circuit model when the electrode resis- 
tances are negligible. 
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FIG. 1. Circuit model for a ferromagnet-insulator- 
-ferromagnet (F-I-F) MTJ device with antiparallel ordered 
moment orientations of the ferromagnets. Each spin chan- 
nel has resistance contributions from the two electrodes and 
from the tunnel junction. In our convention the up-spin chan- 
nel corresponds to minority spins and the down-spin channel 
to majority-spins on the right-hand side of the junction. In 
this configuration the tunnel junction resistances in the two 
spin channels are equal. In panel (a) there is no spin relax- 
ation and up-spin and down-spin channels carry current in 
parallel. In panel (b), spin relaxation is fast enough to main- 
tain equal electrochemical potentials for up and down spins 
throughout the device. In panel (c) equilibrium between up 
and down spins is established when the distance from the 
tunnel junction exceeds the spin-diffusion length of the fer- 
romagnetic electrode, h. Light shading is used for electrode 
resistors in which up and down spins are driven from local 
equilibrium by the inhomogeneous spin-dependent transport 
coefficients of the device. 

For finite electrode resistances and independent con- 
duction in up and down spin channels, we see from the 
circuit model that the local electrochemical potentials 
of up and down spins differ, most strongly at the tunnel 
junction. Whenever this occurs, spin-flip processes in the 
ferromagnet which permit the spin distributions to relax 
to local equilibrium become important. Separate charge 
conservation for up and down spins constrains the parti- 
tioning of current between channels in different parts of 
the device and increases its overall resistance compared 
to the case where spins are free to change their spin-state 
We refer to the resulting increase in resistance as 



the spin-bottleneck resistance of the device. In the fast- 
spin relaxation limit, the spin-bottleneck resistance van- 
ishes and local equilibrium between spins is maintained 
everywhere in the device. In the circuit model this can 
be achieved by shorting the two spin-channels on both 
sides of the junction as illustrated in panel (b) of Figs. |l| 
and |. 
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FIG. 2. Circuit model for a ferromagnet-insulator- 
-ferromagnet (F-I-F) MTJ device with parallel ordered mo- 
ment orientations. In this configuration the two spin channels 
have distinct tunnel junction resistances. The up-spin chan- 
nel corresponds to majority spins and the down-spin channel 
to minority-spins on the both sides of the junction. Light 
shading is used to signal a spin bottleneck as in Fig. |l|. 

The use of a circuit model to represent the transport 
properties of the device is essentially exact in both zero 
and fast relaxation limits, being dependent only on the 
usual locality assumptions. We propose its use in the 
intermediate case where, as we discuss briefly below, lo- 
cal equilibrium between the spin channels is established 
only when the distance from the junction exceeds the 
spin-diffusion length l s . This can be represented in the 
circuit model by shorting the two spin-channels in each 
electrode resistor at the point l s from the tunnel junction 
as illustrated in panel (c) of Figs, [l] and ^. 

The circuit model allows the interplay between spin- 
dependent tunneling and spin-bottleneck effects to be ad- 
dressed in a simple and intuitive manner. Formulas for 
the overall device resistance for both parallel and antipar- 
allel orientations follow immediately from Figs, [l] and || 
using elementary circuit analysis rules. For the interme- 
diate case, < l s < If, we obtain: 
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Note that in a MTJ device the spin-bottleneck appears 
even when the ordered moment orientations are parallel. 
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FIG. 3. Magnetoresistance as a function of the spin dif- 
fusion length calculated using the simple circuit model (full 
lines) and from solutions of linear transport equations for the 
inhomogeneous device (dashed lines). These results are ob- 
tained for a device with comparable electrode and tunnel junc- 
tion contributions to the resistance. The magnetoresistance 
calculated for zero electrode resistance using the Julliere for- 
mula is indicated by dot-dashed lines. 

Typical results for the dependence of overall device 
magnetoresistance on spin diffusion length are shown in 
Fig. ^. These results are for the case of a device with 
comparable electrode and tunnel junction resistances; 
R a pRp/(R a p + R F ) = R T i /2. Spin-bottleneck contri- 
butions to the magnetoresistance decline in quantitative 
importance along with the electrode contribution to the 
overall resistance, but the trends as a function of spin 
diffusion length are similar to those shown here. The 
left panels in Fig. || are for the case t a N a = 2tjNi, while 
the right panels are for the case t a N a — lOijiVj. Ne- 
glecting electrode resistance, these material parameters 
correspond to MR —1/5 and MR = 81/101 respec- 
tively. The curves in Fig. ^| are readily understood using 
the simple circuit model. The dependence of overall re- 



sistance on moment orientation is due completely to the 
central portion of the circuit in which the two spin states 
are not in local equilibrium. For antiparallcl orientations 
this portion consists of two identical resistors in parallel 
with resistance R* = (I s /If)(R f + r f) + r t- For paral- 
lel orientations the two resistors in the central regions are 
not identical and have values R a * = 2(l s /lp)R F + Rp 1 
and RT = 2(l s /l F )R F + R T . Note that the electrode 
contribution to R* is the arithmetic mean of the elec- 
trode contributions to R a * and R % * whereas the tunnel 
junction resistance contribution is the geometric mean 
of the the tunnel junction contributions to R a * and R". 
The property (R*) 2 — R a *R", assuming the Julliere for- 
mula, holds when electrode resistances are neglected and 
guarantees a positive magnetoresistance. This property 
is lost in the general case and MR can be negative. MR 
will tend to be large and positive when either R a * or R" 
is substantially smaller than R*. 

We have considered only cases where the tunneling 
conductance is larger for majority spins, since this ap- 
pears [pUPl to be the case for all systems studied experi- 
mentally. The three panels from top to bottom in Fig. |3| 
are for majority spin resistivities, smaller, equal to and 
larger than minority spin resistivities. All three cases can 
occur [ fl0[ in metallic ferromagnets, depending on the ma- 
terial and the disorder potential. For l s — 0, R A — R p 
equals its Julliere value, and MR is reduced only be- 
cause of the overall resistance increase coming from the 
electrodes. In both top panels, R a */R*, decreases fur- 
ther below 1 as l s increases, leading to an increase in 
MR. For the case R F = R l F illustrated in the middle 
panels, R a */R* increases toward 1 as the electrode con- 
tributions are added and MR eventually changes sign 
as l s approaches lp- As mentioned above MR need not 
be positive in the general case. In the bottom panels, 
R a * /R* increases as l s increases from 0, causing MR to 
change sign, but for larger l s the electrode resistance ef- 
fect begins to dominate, R l * becomes much smaller than 
R* and MR again becomes positive. The minimum MR 
occurs at smaller Z s in the left panel case, because of 
the smaller difference in tunneling resistances. Similar 
simple considerations can be used to obtain an intuitive 
understanding of the magnetoresistance magnitude and 
sign for other material parameter values. 

The electric circuit model, presented above, allows AIR 
values to be understood simply, but is based on an ap- 
proximate representation of the space and spin depen- 
dent local electrochemical potentials in the inhomoge- 
neous MTJ device. To assess its quantitative accuracy 
we have compared its predictions for the overall device 
resistance and magnetoresistance with the results of a 
linear response calculation which assumes only a local 
conductivity for each spin channel Hf7|. In this calcu- 
lation, spin-up and spin-down electrons carry current in 
parallel but spin-relaxation can occur everywhere in the 
ferromagnetic electrodes. Applying Ohm's low locally, 
the current densities for each spin channel are given by 
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where M-m) i s the space and spin dependent electrochem- 
ical Q potential and &-\m = If / '(ARp^). Applying the 
continuity equation to both spin channels and employing 
a relaxation time approximation implies that 
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The non equilibrium magnetization, M(x), is propor- 
tional to the difference between number of spin up 
and down electrons, i.e., M(x) = Mo + Ms Wt^t ( x ) ™ 
N^n^(x)], where jU-m) is the space and spin dependent 
chemical potential, the Fermi energy has been chosen as 
the zero of energy, and Mq is the saturation moment of 
the ferromagnet. M eq (x) is the local quasi-equilibrium 
magnetization towards which M(x) relaxes. If we as- 
sume that the exchange-splitting of the equilibrium fer- 
romagnetic bands Q is not altered by a change in local 
electron density, it follows that M eq (x) = M + /ig(i\^- — 
N±Wl/W + iVj>T + N { /{N^ + N^] and we obtain 
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Eqs. (@) and (§) with the relaxation model (ro) im- 
plies electrochemical potential differences between up 
and down spins which decay exponentially over a length 
scale L inside the electrodes 
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These equations can be solved to determine the partition- 
ing of current between spin channels throughout the de- 
vice. The tunnel barrier enters the calculations through 
the assumption of electrochemical potential drops of 
<57Ij = ARpji and 8ji^ — AR^ji- The other bound- 
ary conditions, needed to solve Eqs.( |^) and (||), express 
the continuity of the current of each spin in the tunneling 
barrier and the recovery of equilibrium spin polarization 
at the ends of the ferromagnetic electrodes |TJ|]. The 
circuit model captures the essence of these equations in 
a simple way. A comparison between the circuit model 
and the macroscopic transport theory outlined above is 
shown in Fig. |[ The two models agree almost exactly 
for l s up to one half of If- The small deviation at larger 
l s results from the fact that circuit model has a cusp at 
l s = If and so is uncapable of describing the smooth 
saturation of the MR curve. The good overall quantita- 
tive agreement between the two models justifies the use 
of the equivalent circuits of Figs. [I] and to understand 



the magnetoresistance of MTJ devices with low overall 
resistances. 
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